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ABSTRACT 

We assess the possibility to detect the warm-hot intergalactic medium (WHIM) in emission and 
to characterize its physical conditions and spatial distribution through spatially resolved X-ray spec- 
troscopy, in the framework of the recently proposed DIOS, EDGE, Xenia, and ORIGIN missions, all 
of which are equipped with microcalorimeter-based detectors. For this purpose we analyze a large 
set of mock emission spectra, extracted from a cosmological hydrodynamical simulation. These mock 
X-ray spectra are searched for emission features showing both the O VII Ka triplet and O VIII Lya 
line, which constitute a typical signature of the warm hot gas. Our analysis shows that 1 Ms long 
exposures and energy resolution of 2.5 eV will allow us to detect about 400 such features per deg 2 
with a significance > 5<r and reveals that these emission systems are typically associated with density 
^100 above the mean. The temperature can be estimated from the line ratio with a precision of ~ 20 
%. The combined effect of contamination from other lines, variation in the level of the continuum, 
and degradation of the energy resolution reduces these estimates. Yet, with an energy resolution of 
7 eV and all these effects taken into account, one still expects about 160 detections per deg 2 . These 
line systems are sufficient to trace the spatial distribution of the line-emitting gas, which constitute 
an additional information, independent from line statistics, to constrain the poorly known cosmic 
chemical enrichment history and the stellar feedback processes. 

Subject headings: cosmology: observations — intergalactic medium — large-scale structure of Universe 
— X-rays: diffuse background 



1. INTRODUCTION 



After more than 10 years, the so-called "missing 
baryon" probl em is far from being solved . The orig- 
inal claim of iFukugita. Hogan. fc Peebles! (|1998| ) that 
the observational census of cosmic baryons in the lo- 
cal Universe falls short of the expected cosmologi- 



cal mean (IBurles fc Tvtlerl [19971: IBarlow fc Tvtleil HIM 
iSpergel et al.l2003UKomatsu et al.ll2009ri has gaine d even 
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more statistical signifi cance ( Fukugita fc Peebles! 120041 : 
lDanforth~fc Shull 2005) and it is now generally accepted 
that a large fraction of the baryons (~40 %) is not ac- 
counted for by observations at z ~ 0. This fact is 
even more surprising if one considers the fact that the 
baryon density in the Ly-q forest at redshift z =2-3 
(jWeinberg et al.l Il997t iRauchl 11998ft accounts for about 
all expected baryons whereas the nucleons in the Ly-a 
forest at z ~ contributes to less than 30 % of the cos- 
mic mean (|Bregmaril2007i) . 

Although baryons provide a small contribution to the 
cosmological mass-energy density, they are key players in 
forming stars and galaxies and determining the cosmo- 
chemical evolution of the diffuse gas. In addition, being 
the only component that interacts with the electromag- 
netic radiation, they are the sole probe to investigate the 
physics of the formation and evolution of visible cosmic 
structures and infer the co-evolution of the underlying 
dark matter distribution. 

Cosmological numerical simulations have provided a 
solution to this apparent paradox suggesting that miss- 
ing baryons at z < 1 are distributed in a network of fil- 
aments with temperature 10 5 -T0 7 K with density signif- 
icantly above the cosmic mean but below that typically 
found in virialized structures like galaxy clusters. This 
phase is called warm- hot intergalactic medium (WHIM) . 
These results, which constitute one of the rare examples 
in cosmology in which theory leads observations, explain 
why these baryons have escaped observations so far and 
indicate the best strategy to detect them. Indeed, the 
WHIM is expected to be so highly ionized that it can 
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only be seen in the UV and X-ray bands but the signal 
should be very weak because of its relatively low density. 
Therefore it is not surprising that current observations 
in the UV and X-ray bands provide weak constraints to 
the WHIM models. 

In the UV band, O VI absorption lines (A = 
1032, 1038 A) have bee n detected along the lines of 
sight to ~ 50 AG N (e.g., IDanforth fc Shulll I2005L [20081: 
iTripp et all I2008T ) . Their cumulative number per unit 
rcdshift as a function of the line equivalent width (here- 
after the dN/dz statistics) are consis tent with the pre- 
dictions of numerical models (see e.g.. lCen fc Fand[2006t 
iBranchini et al.ll2009D . However, it is debated how reli- 
ably the observed O VI absorption lines trace the WHIM 
since their width seems too small to be produced by 
warm-hot material. This means that either the observed 
lines trace the local Ly-a forest rather than the WHIM or 
that the gas is two-phased: a photo-ionized gas respon- 
sible for the Ly-a absorption and shock-heated WHIM 
responsible for the O VI lines. Recent numerical sim- 
ulations carried out to investigate the physical proper- 
ties of the IGM have helped to clarify the issue (????). 
These simulations have shown that the thermal state of 
the IGM at z ~ is significantly affected by a num- 
ber of physical processes like galactic wind, AGN feed- 
back, metal line cooling, photoionization by background 
radiation, sub-resolution turbulence. Although it is non 
trivial to compare the results of different works and not 
all these works seem to converge to the same results, 
it is remarkable that two of them (??) indicate that 
about a third of the OVI absorbers are photo-ionized 
and trace the Ly-a forest at z ~ while the rest of them 
trace the WHIM. An upper limit can nevertheless been 
worked out: if all gas is multi-phased, i.e., in the limit 
in which all the observed O VI lines are produced by the 
WHIM, then the O VI can probe no more than 7-10 % of 
the baryons (IDanforth fc Shulll [20051: ITripp et al.l [2^51 
iDanfortbl 12009^ With the inclusion of the thermally 
broadened Lyq absor bers ( BLAs) recently o b served by 
IDanforth et al.l (|2010f ) and iNaravanan et al.l (|2010al lbl). 
these estimates are boosted up by a factor of ~2, i.e., UV 
observations are capa ble of probing ~20% of the baryons 
in the local universe (IDanforth et al.|[2010[ ). Neverthe- 
less, whereas O VI might be a useful signpost for a non 
negligible fraction of the missing baryons, the bulk of the 
WHIM cannot be investigated through UV spectroscopy. 
X-ray band observations are required. 

In the soft X-ray band, the residual unresolved back- 
ground after point-sources subtraction possibly consti- 
tutes the second observational evidence of the WHIM. 
So far the tightest c o nstrain ts have been obtained by 
iHickox fc Markevitchl (j2007al lbh after removing all con- 
tributions associated to X-ray, optical, and IR sources 
in the Chandra Deep Fields. Their analysis has shown 
that the residual background in the 0.65-1 keV energy 
band cannot be accounted for by AGN below the resolu- 
tion limit. Instead, it appears to be consistent with the 
expected integrated emission contributed by the WHIM, 
as predicted by models based on h ydrodynamical simu- 
lations (e.g.. lRoncarelli e t al. 200j|). The intensity of the 
diffuse X-ray background constitutes a weak constraint 
on the WHIM properties. Additional constraints can be 
obtained by identifying the WHIM signature in the an- 



gular correlation propert ies of the X-ray backg round, like 
in the recent analysis of iGaleazzi et al.l (|2009l ). 

Hydrodynamical simulations show that the best way 
to investigate the WHIM is by observing the X-ray emis- 
sion or absorption lines of highly ionized element like 
C, N, O, Ne, Mg, and Fe. Deep X-ray spectroscopy 
has indeed been performed to detect such lines. i,From 
a theoretical viewpoint the possibility of detecting and 
studying the WHIM in the absorption spectra of X- 
ray bright objects like AGN or Gamma-ray bursts has 
been thoroug hly discussed both f r om a general per- 
spective (e.g.. iKraytsov et all l2002t iKlvpin et al l 120031 : 
IChen et all 12003 iViel et"aiT[2005T ICen fc Fand 12009) 



and by considering o bservations carri e d out with some 
speci fi c ins t rument (IViel et al.l [20031: lY oshikawa et alj 
2001 [2001 iKawahara et al.1 12009 iPaerels et al.1 120081 : 
Branchini et al J [20091). ; From the observational side, 



the analysis of the absorption spectra has lead to sev- 
eral claimed WHIM detections (e.g.. IFang et al.1 120021: 
iMathur et all l2003t iNicastro etall I2005D whose statis- 
tical si gnificance has been questioned by subseq uent 
studies (jRasmussen et al.1l2007t iKaastra et alJ[2(509) . re- 
flecting the fact that these measurements are at the 
limits of current instrumental capabilities. The at- 
tempts to detect absorption lines associated to dense 
parts of large-scale structure lead s to marginal detec 
tions with less than 4cr significance (iFuiimoto et all 



[ Takei et al.1 120071: iBuote et al.l 120091 : IFang et al.1 
iZappacosta et al* 120101 ) 



2004; 



2010; 



Studying the WHIM in absorption offers advantages 
and disadvantages. The advantage is that line detection 
depends on the source flux with no need to worry about 
foreground contributions. In addition, the line strength 
depends on the density of the intervening absorber, al- 
lowing one to probe regions of moderate overdensity. The 
obvious drawback is that one can only probe the WHIM 
along the few lines of sight to X-ray bright objects, mak- 
ing it very difficult to trace its spatial distribution (but 
see IViel et ail 120031 for a possible way of probing the 
WHIM clustering in absorption). 

On the other hand, studying the WHIM in emis- 
sion is considerably harder but potentially more re- 
warding. The main difficulty is represented by the 
fact that line detection depends now on the capabil- 
ity of modeling the various contributions to the soft X- 
ray foreground, which, unfort unately, are poorly con- 



strained by observational data ( McCammon et al.| | 2002 t 
Smith et all 120051: iHenlev et al.N2007l: iSmith et aLlboOTt 



Henlev fc S hclton 2008; Gupta et al.ll2009HYoshino et atl 



20091) . It is therefore not surprising if only a few 
WHIM d etections have been claimed in continuum 
emission (IZappacosta et "ahl 120051: iMannucci et al. l 120071 : 
iTakei et al.1 120071 IWerner et al.l l2008f h The main re- 
ward is constituted by the po ssibility of tracing the 
3D distribution of the WHIM (jCen fc Fand 12009) and 
quantifying its spatial and angular correlation properties 
(jUrsino fc Galeazzi|[200l lUrsino et al.ll2010bD . In addi- 
tion, since X-ray emissivity scales with the square of the 
gas density, we expect the emission signal to be domi- 
nated by high density regions, effectively complementing 
the absorption studies. 

In this work we will investigate the possibility of detect- 
ing, studying and characterizing the WHIM in emission. 
In order to observe the weak WHIM emission features a 
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detector with a large grasp (effective area x field of view) 
is mandatory. In addition, to separate the WHIM lines 
from other spectral features contributed by other sources 
a good energy resolution is also required. A number of 
recently-proposed X -ray satelli te missions such as DIOS 
• Qhashi et all l2010h. EDGE (|Piro et al.l [2009h . Xenia 



Burrows et al.ll20Toh , and ORIGIN meet these require- 



ments. We will assume observations with the detectors 
on board the proposed EDGE and Xenia missions. The 
detailed char acteristics o f thes e two X-ra y satellites are 
presen ted in IPiro et all (|2009lf TS L and in iBurrows et alJ 
(|2010T l . For our purposes, the most relevant instrument 
is the cryogenic microcalorimeter-bascd spectrometer ar- 
ray: CMS (Cryogenic Imaging Spectrometer). CRIS has 
the energy range of 0.2-2.2 keV (goal 0.1-3 keV), the en- 
ergy resolution of Ai?=2.5 eV (goal 1 eV) at 5.9 keV, 
the FOV of 0.9° x 0.9° (goal 1.1° x 1.1°), the effective 
area A = 1000 cm 2 (goal 1300 cm 2 ) at 0.5 keV, and the 
angular resolution (HPD) of A6 = 4' (goal 2.5'). DIOS 
spectrometer has a similar energy resolution but a factor 
of 10 smaller grasp than CRIS onboard EDGE or Xenia. 
ORIGIN is planned to have a better angular resolution 
< 30", at the expense of energy resolution (AE=3— 5 
eV) . In this work we refer to the adopted baseline of the 
EDGE and Xenia spectrometer as a reference case and 
explore the effect of a different energy resolution, chosen 
in the range [1-7] eV. In contrast, we do not explore the 
angular resolution better than 1.3', because the typical 
angular size of the WHIM above detection threshold is 
> 2'-3'. 

The layout of this paper is as follows. In §[2] we intro- 
duce the WHIM model used in this work from which we 
obtain the simulated surface brightness maps and emis- 
sion spectra analyzed in this work. In §[3] the problem of 
detecting the WHIM emission lines and computing their 
relevant statistics is addressed. In §[4] we study the phys- 
ical properties of the regions in which the WHIM can be 
detected in emission and assess the possibility of measur- 
ing the gas temperature from line ratios. In §[5] we inves- 
tigate the possibility of tracing the spatial distribution 
of the WHIM. In § [SI we quantify the impact of energy 
resolution and line contamination form different sources. 
In § [3 we address the question of to what extent we can 
discriminate among different models of stellar feedback 
and metal diffusion, using emission studies. Finally we 
discuss the main results and conclude in § [5] 



2. MODELING THE WHIM 



Hydrodynamical simulations provide the 
of studying the prop ert ies of the WHIM (see, 
Cen fc Ostriken Il999allb|: IDave et alJ [200lt IChen et alJ 



20031: IBorgani et al 



best way 
e^ 



20041: 

Qppenheimer & Pavel 120081 iBertone et alJ 12010 : 
Tornatore et al.ll2010( ). Their analysis has revealed that 
the missing baryon problem finds its natural solution in 
a standard ACDM framework in which the WHIM is 
heated up to temperature of 10 5 -10 7 K mostly by means 
of hydrodynamical shocks resulting from the build-up 
of cosmic structures at scales that are entering the 
nonlinear regime of density fluctuations growth. Addi- 
tional heating/cooling mechanisms like stellar and AGN 
feedback, galactic superwinds, and radiative cooling, 

15 Also sec http://projects.iasf-roma.inaf.it/EDGE 



iCen fe Ostrikeif 
IBertone et alJ 



2006 



and ionization by background radiation further modify 
the thermal state of the gas. Desp ite the recent advances 
in modeling these mechanisms ( Cen fc Ostrikerl 120061: 
Cen fc Fangj 120061: IBertone et al.l [20101 : iTornatore et al l 



20101: iShen et alJI^OlOb ?: ?: ?: ?) we still do not have 



a firm physical understanding of the effects of stellar 
feedback, chief among which is metal enrichment. 

Indeed, the metal enrichment is so important in deter- 
mining t he observational p ropert ies of the WHIM that, 
following IBranchini et al.l (|2009l ). we prefer to specify 
gas metallicity after the simulation run (i.e in the post- 
processing phase) rather than self-consistently during 
the numerical simulation. The resulting WHIM model 
has been shown to fulfill all existing observational con- 
straints ranging from the dN/dz of the O VI and O VII 
absorption lines to the surface brightness of the diffuse 
X-ray background. So, although effective, the WHIM 
model used in this work is not self-consistent. However, 
given the existing uncertainties in predicting metal abun- 
dances, in this work we will not rely on a single, sophis- 
ticated WHIM model. Instead, in an attempt to account 
for theoretical uncertainties, we will consider two WHIM 
models that have been obtained from the same hydrody- 
namical simulations but adopting different prescriptions 
for the gas metallicity. 

Our WHIM models have been obtained f r om th e hy- 
drodynamical simulation of IBorgani et alJ (120041) per- 
formed with the GADGET-2 Lagrangian code (jSpringell 
120051) in a computational box of size 192/j" 1 comov- 
ing Mpc, loaded with 480 3 Dark Matter and 480 3 gas 
particles. Dark matter and gas particles have masses 
of 4.62 x 10 9 h~ 1 M Q and 6.93 x lO 8 /^ 1 M Q , respec- 
tively. The Plummer-equivalcnt gravitational softening 
is e=7.5 hr 1 kpc at z = 0, fixed in physical units be- 
tween z = 2 and z = 0, i.e., the redshift interval relevant 
for the study of the WHIM. The background cosmolog- 
ical model is a flat ACDM with il A = 0.7, fi b = 0.04, 
H = 70 km s" 1 Mpc 
normalization a % = O.i 
simple recipes to account for the main non-gravitational 
heating and cooling mechanisms, namely: i) the star 
formation process that is treated by adopting a sub- 
resolution multiphase model for the interstellar medium 
(|Springel fc Hernquistl 120031 ). ii) the feedback from SNe 
that includes the effect of weak galactic outflows and Hi ) 
the radiative cooling assuming zero metallicity and heat- 
ing/cooling from t he ph oto-ionising UV background by 
lHaardt fc Mad au (1996). Finally, to construct our mock 
emission spectra we have considered all available simula- 
tion outputs out to z = 0.5. 

As mentioned above, the gas metallicity has been spec- 
ified in the post-processing phase and, to account for the 
scatter in model predictions, we have explored two differ- 
ent scenarios. In the first one, which we call model Bl, 
we have adopted the deterministic metallicity-density 
relation Z = min ( 3, 0.005(1 + 8„) 1/ 2 ), where S g is 
j°gas/(/3gas) — 1, that iCroft et all (|2001l ) have proposed 
to match the average relation measured in numerical ex- 
periments. This is a conservative model that typically 
underestimates the line emissivity since it ignores the 
large scatter in the Z — S g relation which would enhance 
the contribution of the line emission from regions charac- 
terized by moderate gas overdensity, like the WHIM. For 



1 (h = 0.7), and power spectrum 
This numerical experiment uses 
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this reason, the Bl model should be regarded as a gen- 
erous lower bound to our detectability estimates. In the 
second model, dubbed B2, we do account for the large 
scatter in the metallicity-density relation. This scatter 
reflects the fact that the metal abundance at a given loca- 
tion is not simply determined by the underlying gas den- 
sity but also by the thermal and chemical history of the 
gas particles. Since it would be hard to model this scat- 
ter analytically we resort to hydrodynamical simulations. 
For this reas on we have enforced the s ame Z — S g relation 
measured by iCen fc Ostrikerl (|1999al fbl) in their numeri- 
cal experiment. In practice we use the de nsity-metallicity 
scatter plot measured at z = in the ICen fc Ostrikerl 
(|1999bft simulation to derive a 2D probability distribu- 
tion which is then fed into a Monte Carlo procedure to 
assign metallicity to gas p articles with known density 
in the IBorgani et all ( 20041 ) simulation. The result is a 



metallicity model with the same Z — S g relation plotted 
in Fig. 2 of ICen fc Ostrikerl (|1999aD and characterized 
by a 1 a scatter of ~ 0. 3 dex. Note that we did not 
use the metallicity of the IBorgani et al.l (|2004l ) simula- 
tion itself since the metal diffusion mechanism built-in 
the SPH code reproduces the correct metal content in 
high density regions, like the intracluster medium, but 
systematically underestimates the metal abundance in 
the typical WHIM environments. 

We stress the fact that the metallicity assigned in the 
post-processing is not self-consistent with the thermal 
state of the gas which, in the simulation considered in 
this work, is determined ignoring the role of metal cool- 
ing. In fact, recent numerical experiments show that the 
thermal state of the diffuse gas can be significantly al- 
tered by metal cooling (????). However, including the 
effect of metal cooling is beyond the scope of this paper. 
Instead, we attempt to bracket the uncertainties that 
derives from our simplified by adopting the two very dif- 
ferent density-metallicity relations, being confident that 
the deterministic density-metallicity relation will provide 
very conservative predictions for the detectability of the 
WHIM. 

To determine the ionization balance we have assumed 
pure collisional ionization equilibrium (CIE) and not hy- 
brid collisional and p hotoionization equilibrium as in 
IBranchini et~aTI (|2009f ). The use of the CIE assumption 
is justified by the fact that in this work we are inter- 
ested in emission spectra whose intensity scales with the 
squared density of the gas and thus it is heavily weighted 
towards high density peaks wher e the ionisation i s basi - 
cally determined by collisions. IBranchini et al.l ([20091 ) 
were interested in absorption lines which typically trace 
the WHIM in regions of moderate overdensity where 
photo-ionization cannot be neglected. Here we focus 
on WHIM lines For this reason we have assumed CIE 
and used the APEC therma l model included in the code 
XSP EC dSmith et al.l I2001D to compute emission spec- 
tra (|Ursino fc GaleazzT 20061) . The resulting mock X- 
ray spectra contains all metal lines but here we focus on 
the strongest ones only: the O VIII Lya line (E =0.653 
keV) and the O VII Ka triplet (E= 0.561, 0.569, and 
0.574 keV). 

Model B2 agrees with all available observational 
constraints and, from a theoretical viewpoint, pre- 
dicts a phase-space diagram (shown in Figure 6 of 
IBranchini et al.l (|2009l) ) which is remarkably close to that 




0.65 
Energy (keV) 

Fig. 1. — Emission spectrum in a 2.6' X 2.6' area taken in a 1 
Ms observation with Xenia CRIS, assuming an energy resolution 
AE = 1 eV. Black: Sum of the Galactic Foreground and unre- 
solved extragalactic background. The Galactic O VII triplet and 
the O VIII Kq line at z = are flagged. Red: Contribution from 
the extragalactic gas. The O VII and O VIII lines at z = 0.033 arc 
indicated in the plot. 



of ICen fc Ostrikerl (I2006D (shown in Figure 6 of that 
paper), a fact that further increases our confidence on 
model B2. Therefore in this work we consider B2 as the 
reference model and regard predictions from model Bl 
as a generous lower bound. 

One final caveat for model B2: this model is known to 
overestimate the X-ray emissivity in high (<5 ff > 10 3 ) den- 
sity envi ronments, i.e. emiss ion from gas in groups and 
clusters (jUrsino et al.ll2010al ). While the precise reason 
for this systematics is unknown, a plausible explanation 
is represented by the fact that in our model we have as- 
sumed primordial composition in the cooling function, in- 
stead of updating the cooling rate accordin g to the evolv- 
ing m etallicity of the gas. The analysis of iBertone et al.1 
(|2010( ) suggests that this would artificially increase the 
emissivity in regions with high density. Brighter emission 
from group- like environments increases the chance super- 
position from strong emission lines, like those of the Fe 
L complex, that may contaminate or even outshine the 
weaker O lines, artificially reducing our theoretical esti- 
mate of WHIM detection. We will quantify these effects 
in § 16.11 We again stress the fact that the inclusion of 
metal line cooling might have an impact on both the line 
surface brightness and the thermal state of the gas. We 
qualitatively discuss the expected impact of metal cool- 
ing in § [H 

2.1. Mock emission spectra 

The detailed procedure use d to produce the mock spec- 
tra is described in detail in lUrsino et al.l (|2010bf ) , here 
we provide a brief summary. To compute the mock X- 
ray spectra we have considered all gas particles within 
a light-cone extracted from the hydrodynamical simu- 
lation, out to z = 0.51. The light cone is formed by 
stacking the simulation outputs corresponding to 7 red- 
shift intervals, each one of them corresponding to a co- 
moving depth of 192/i _1 Mpc. We do not consider the 
WHIM at z > 0.5 since in that case the O VII lines would 
be redshifted to energies at which a number of emission 
lines due to L and M transitions of heavier metals make 
identification of O lines difficult. In the stacking pro- 
cess, the simulation cubes that constitutes the individual 
rcdshift intervals were randomly shifted and rotated to 
avoid periodic replicas of the same large scale structures. 
Then, after selecting a random observer within the cube 
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at z — we generate a single cone of view of 5.5° x 5.5° 
and tag gas particles with their redshift, angular position, 
density, and temperature. Particles with temperature in 
the range 10 5 — 10 7 K and overdensity S < 1000 are 
classified as "WHIM phase" . Among the remaining gas 
particles, those that provide a significant contribution to 
line emission are either in a "hot phase" (S > 1000 and 
T > 10 7 K, roughly corresponding to gas within clusters 
of galaxies) or in a "dense phase" (S > 1000, T < 10 7 
K, typically associated to galaxy groups). Their con- 
tribution to the X-ray emission is included in our mock 
spectra. The contribution from particles colder than 10 5 
K is negligible. In the mock spectra we did not include 
the contribution from gas at z > 0.51. Emission lines 
(like the Fe L complex in galaxy groups) from high-z 
gas potentially contaminates spectra in the energy range 
in which we search for WHIM lines. We estimate the 
impact of contamination by spurious lines in § 16.11 

To account for the finite angular resolution of the in- 
strument, the FOV of the light cone is divided into 
256 x 256 pixels of 1.3' x 1.3', i.e. smaller than the 
planned CMS resolution but possibly consistent with 
alternative instrument designs. The mock X-ray spec- 
trum within each angular resolution element was ob- 
tained by adding together the individual spectra from 
each gas particle within the pixel plus that of gas par- 
ticles in the contiguous pixels, weighted according to 
the smoothing kernel. The mock spectra were then 
sampled with an energy resolution of 1 eV to mimic 
the energy resolution of the instrument. The effect of 
the Gala ctic absorption has been inclu ded using the 
model of iMorrison fe McCammonl (|1983|) in which we 
have adopted typica l high latitude column de nsity value 
of 1.8 x 10 cm - 2 (jMcCammon et all 120021) . We have 
assumed the same meta l solar abundance as in the 
lAnders fe Grevessi (|1989f ) models. The final result is 
a suite of 256 x 256 x 7 mock spectra that we search 
for detectable WHIM lines. A typical example of mock 
spectrum is shown in Fig. [T] The red curve represents 
the contribution by extragalactic gas in a light cone with 
FOV of 2.6' x 2.6' (i.e. the sum of 2 x 2 x 7 individual 
mock spectra). Photon counts have been computed in 
AE = 1 eV energy bins assuming a 1 Ms observation 
with CRIS. The spectrum is convolved with the detec- 
tor response matrix of 1 eV energy resolution (FWHM) . 
Emission lines from different metals are clearly visible. 
In particular we have flagged those lines corresponding 
to a line system at z = 0.033. The O VII triplet and the 
O VIII recombination line are clearly visible. The black 
curve shows the sum of the contributions of the Galactic 
foreground an d the unresolved ex t ragala ctic background 
modeled after iMcCammon et all ()2002f ) . The Galactic 
O VII and O VIII lines at z = are clearly seen in corre- 
spondence of their rest frame energy. 

3. ANALYSIS OF THE MOCK SPECTRA 

To assess the possibility of detecting the WHIM in 
emission and constraining its thermal state with next- 
generation instrument we consider spectra in larger pix- 
els (2.6' x 2.6') similar to the CRIS angular resolution 
(goal). For this purpose we search for emission lines in 
each of the 128 x 128 x 7 mock spectra using a two step 
procedure. (1) We identify the O VII (triplet) and O VIII 
lines by searching for local maxima in the appropriate 



energy range of the spectrum. (2) we compute the sur- 
face brightness of each line by summing the flux over 
energy bins redward and blucward of the line maximum 
and stop whenever the surface brightness SB t /j in the 
bin drops below 7 x 10~ 3 ph s _1 cm~ 2 sr _1 which is 
smaller than the typical line detection threshold for a 1 
Ms observation with CRIS, as we will show § 13.21 A sim- 
ple example serves to clarify the procedure. To identify 
O VIII lines produced in the redshift range z = [0, 0.065] 
we search for local maxima in the energy range [0.613, 
0.653] eV, where the lower bound is the redshifted energy 
of an O VIII photon emitted at z = 0.065 The ccntroid 
of the line is identified with the maximum and its inten- 
sity is obtained by summing over all contiguous energy 
bins with surface brightness above that of the selected 
threshold. This strategy, which minimizes the chance 
contamination from gas at different redshifts, cannot be 
applied to real spectra since we lack the information of 
the redshift of the gas responsible for the line emission. 
This effect is discussed in § 16.11 

Detectable WHIM lines are rare and it is unlikely 
to find two of them in the same mock spectrum. On 
the other hand, when approaching the surface bright- 
ness threshold, several weak emission lines start to ap- 
pear. They are typically clustered around stronger lines 
and thus are likely to be physically associated to the 
same line emitting element. Considering all these lines 
would spoil line statistics since they would oversample 
the same line-emitting regions. To circumvent this prob- 
lem we have adopted a simple criterion: we merge to- 
gether all lines that are separated by a number of en- 
ergy bins smaller than some minimum amount ribin- In 
our analysis we have adopted Ubm = 8 (corresponding 
8 eV) after having checked that the number of lines with 
0.1 ph s _1 cm~ 2 sr _1 (strong enough to be hopefully de- 
tected by next-generation instruments) depends neither 
on nbin nor on SB t ^. 

To confirm that our B2 model fulfill the current con- 
straints on the soft X-ray diffuse background, we have 
compared the total surface brightness predicted by the 
model in the [0.65, 1] keV ba nd with that measured by 
iHickox fe Markevitchl (|2007al) after removing the contri- 
bution from all known sources. According to the model, 
the WHIM contribution to the surface brightness con- 
tributed is (3.6±0.3) x 10~ 13 erg cm~ 2 s" 1 deg~\ safely 
below the observed value of (1.0 ± 0.2) x 10~ 12 erg cm~ 2 
s" 1 dcg -1 . The more conservative Bl model, which pre- 
dicts a lower surface brightness is also obviously in agree- 
ment with observations. 

3.1. Line Statistics 

After having identified all O VII and O VIII emission 
lines in the mock spectra we can compute their cumu- 
lative number per unit redshift as a function of the line 
surface brightness. The top panel of Figure [5] shows the 
case of the O VII triplet whose surface brightness is the 
sum of the three lines. The black continuous curves refer 
to model B2 while the red dashed curves refer to model 
Bl. For each model, the upper curve accounts for the 
gas emission from all gas particles along the line of sight 
and the lower curve considers only the contribution from 
the "WHIM phase" (see §EJ). In the B2 model the lines 
contributed by the warm hot gas constitutes a significant 
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TABLE 1 
Emission line detections 



/lino dNovu/dz dNoviu/dz dNovu+OVlll/dz TVovii+ovm per deg 2 

B2 1.0 0.07 4.2 (3.1) 2.6 (1.9) 2.4 (1.6) 639 (426) 

B2 0.1 0.48 2.1 (0.6) 1.5 (0.5) 1.1 (0.2) 293 (53) 

Bl 1.0 0.07 3.3 (1.6) 2.1 (1.0) 1.8 (0.7) 479 (186) 

Bl 0.1 0.48 1.4 (0.03) 1.1 (0.05) 0.7 (< 10" 3 ) 186 (0) 



Note. — Column 1: WHIM model. Column 2: Exposure time (Ms). Column 3: Minimum 
line surface brightness required for a 5<r detection (ph s -1 cm -2 sr -1 ). Column 4: Expected 
number of O VII detections per resolution element and unit redshift Column 5: Expected number 
of O VIII detections per resolution element and unit redshift. Column 6: Expected number of 
simultaneous O VIII and O VIII detections per resolution element and unit redshift. Column 
7: Expected number of simultaneous O VIII and O VIII detections per square degree due to 
gas within z = 0.5. All estimates assume an angular resolution of 2.6' X 2.6'. The numbers in 
parenthesis indicate lines contributed by the WHIM. 




Surface brightness (ph s~ 1 cm- 2 sr 1 ) 

Fig. 2. — Cumulative number of emission lines per unit redshift 
as a function of line surface brightness. Each panel refers to O VII 
triplet (top), O VIII Lyo (middle) or both O VII triplet and O VIII 
Lya (bottom). Black, solid curves indicate model B2, while red, 
dashed curves model Bl. For each model, upper curve shows all 
the gas and lower curve refers to "WHIM phase" gas only. 

fraction of the total. However, this fraction decreases 
with the surface brightness of the lines. It ranges from 
~ 100 % for the weakest line to ~ 5% for the brightest 
one. This result shows that the possibility of detecting 
the WHIM emission lines increases when decreasing the 
line detection threshold and clearly illustrates the need 
of increase the sensitivity of next-generation instruments. 
Model Bl predicts ~ 30% less lines than model B2, as 
expected. The difference among the models is almost in- 
dependent of the line surface brightness. For a reference 
line surface brightness of 0.1 ph s _1 cm -2 sr -1 (match- 



ing the typical detection threshold for next-generation 
X-ray spectrometers as we will show) the expected num- 
ber of detected lines is of order of unity, indicating that 
the chance of multiple line detections in a pixel is rather 
small. 

The middle panel of Figure [2] shows the cumulative 
dN/dz statistics for the O VIII line. All considerations 
made for the O VII lines apply to the O VIII case ex- 
cept for the fact that, for a given surface brightness the 
expected number of O VII triplets is larger than that of 
O VIII lines. This difference becomes smaller if one con- 
siders the brightest O VII lines rather than the triplet. 
The bottom panel shows the cumulative dN/dz of those 
line systems that exhibit both the O VII triplet and 
O VIII line as a function of the total surface brightness. 
The curves are similar to those of the O VIII lines (mid- 
dle panel) except that at the bright end. This is due to 
the fact that presence of an O VIII line almost invari- 
ably guarantees that of an O VII triplet with similar (or 
larger) surface brightness. 

Our resu lts can be compared with the theoretical pre- 
dictions of iCen fc Fand (|2006l ) who used more sophisti- 
cated models that account for departures from ioniza- 
tion equilibrium and the effects of galactic supcrwinds. 
Their predicted cumulative number of emission lines (see 
Fig. 12 of the paper) is similar to that predicted by model 
B2. In particular the agreement is remarkable for the 
O VIII lines, while we seem to over predict that of O VII 
triplets. In addition, we note that the difference between 
model Bl and B2 is smaller than the scatter in their 
model predictions, indicating that we are being generous 
in our attempt of bracketing theoretical uncertainties. 

3.2. Expected WHIM detections 

To detect O VII and O VIII emission lines one needs 
to discriminate the genuine, weak emission signal from 
the astrophysical foreground and background which rep- 
resents the main source of noise, typically exceeding the 
instrumental background. The problem of detecting the 
WHIM emission lines is analogous to that of detecting a 
low-surface brightness source with a typical angular ex- 
tension is a few arcmin, i.e., the typical angular size of 
a WHIM filament at z ~ 0.3, The signal to noise is then 
optimized by measuring spectra with an angular resolu- 
tion comparable to that of the sources over large FOV 
to detect as many WHIM lines as possible. In practice, 
one wants to perform spatially resolved spectroscopy over 
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large areas with an angular resolution of few arcmin. The 
CRIS detector proposed for the Xenia mission and that 
we are considering as a prototype for next-generation 
spectrographs fulfills these requirements. 

To investigate the possibility of detecting the WHIM 
through line emissions and assess the statistical signifi- 
cance of the detections one needs to quantify the noise 
which is mainly contributed by two sources: the cosmic 
X-ray background (CXB, hereafter) produced by unre- 
solved extragalactic point sources, and the foreground 
emission of our Galaxy. 

The CXB is contributed by all X-ray sources (mostly 
AGN) below the source detection threshold. Its inten- 
sity depends therefore on the characteristics of the de- 
tector and on the exposure time. A conservative upper 
limit to the CXB can be obtained fr om the XQC rocket 
experiment (jMcCammon et al.l[200l ) since this estimate 
includes the contribution of all point sources in the FOV 
of ~1 sr. In this case the CXB surface brightness is 
/cxb ~ 30 ph s _1 cm~ 2 sr _1 kcV -1 . This value can be 
reduced by removing all resolved sources. One can re- 
move sources previously detected in some other surveys. 
However, efficient removal would be possible only in ar- 
eas covered by deep surveys, which are typically small. 
The alternative is to perform simultaneous deep imag- 
ing and resolved spectroscopy of the same areas, a strat- 
egy that would be possible if both a CCD-like imager 
and a 2D spectrograph with similar FOVs were avail- 
able. This is precisely the case of of the Xenia mission 
concept which should be capable of resolving and remov- 
ing ~ 1000 AGN per square degree down to a limiting 
flux of 2.5 x 10~ 16 erg s" 1 cm" 2 in the 0.5 - 2 kcV band. 
If AGNs were isotropically distributed, then a spectro- 
graph with ~ 4' angular resolution would have all its pix- 
els contaminated by AGN emission. Taking into account 
the a ngular clustering of the AGNs (Vikhlinin & For man 
|1995[) . this contamination would affect ^63% of the pix- 
els if the angular resolution of the instrument is improved 
to 2.5' (i.e., the CRIS goal). It is clear that the best 
compromise between low AGN contamination and large 
number of available pixels depends on the actual charac- 
teristics of the detector and on the availability of a CCD 
imager. For the sake of simplicity in this work we adopt 
the conservative solution of keeping all pixels and thus 
we assume /cxb = 30 ph s _1 cm -2 sr _1 keV -1 . 

The Galactic foreground is basically contributed by 
three sources: the so-called local bubble, the Galactic 
halo, a nd the solar wind charge exchange (SWCX). Fol- 
lowing iMcCammon et all (|2002T ) we describe the fore- 
ground emission as the sum of absorbed thermal model 
(wabs and apec model in the XSPEC package respec- 
tively) and unabsorbed thermal model (apec). The for- 
mer component represents the Galactic halo, while the 
latter the sum of local bubble and SWCX. One may think 
that this model is too crude, since the SWCX contribu- 
tion, which is non thermal, is expected to outshine that 
of the local bubble. Moreover, the variable nature of the 
SWC X and the lack of strong observational cons traints 
(e.g. JHenlev k Sheltonll2008l : lFuiimoto et al.ll2007t ) make 
it difficult to model its contribution. However, for the 
purpose of this work, this approximation is reasonable 
for two reasons. First, it has been shown that the sum of 
two thermal components, although inaccurate, is quite 
effective in reproducing observations, at least within the 
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Surface brightness (photons s' 1 cm 2 sr 1 keV" 1 ) 

Fig. 3. — Fraction of energy bins (each bin with 1 eV width) in 
which the measured surface brightness is above the value indicated 
on the X-axis. Different curves represent different energy resolu- 
tion (FWHM) of the instrument: 1 eV (red dashed), 3 eV (green 
dotted), 5 eV (blue dot-dashed), and 7 eV (magenta dash dotted). 
The black, solid curve shows the case of no convolution with the 
detector response. 

energy resolution of CCD detectors (jGaleazzi et alJ|2"o"07t 
lYoshino et al.ll2009T ). Second, in the analysis of the real 
spectra one ignores the energy ranges in which Galactic 
emission lines are present. These lines include Ly/3, Ly7, 
and LyJ transitions from H-likc C and N atoms, which 
should constitute strong features o f the SWCX emis- 
sion. I n fac t, all SWCX lines liste d bv lKoutroumpa et al.l 
(|2006| ) and lBodewits et al.1 (|2007t ) with a contribution of 
>a few% of O VII triplets are in energy intervals that we 
ignore when we search for the WHIM lines in § 16.11 

Since the energy spectrum of the Galactic foreground 
is characterized by prominent emission lines (see Fig. 2) 
and we are interested in identifying those produced by 
the WHIM, Galactic foreground flux cannot be estimated 
by the averaged flux over the wide energy range (e.g., 
0.4-0.7 keV). More realistically, we would need to restrict 
the search to those energy bins that are not heavily con- 
taminated by Galactic lines. Once again, when adopting 
this strategy one should compromise between the need of 
keeping as many energy resolution elements as possible 
and that of minimizing the Galactic signal. The solu- 
tion can be found by looking at Figure [3] that shows the 
fraction of energy bins in which the measured foreground 
surface brightness is above the value indicated on the X- 
axis. The shape of the curves that refer to the case of 
no convolution with the detector response (black, solid) 
and of 1 eV energy resolution (red, dashed) is charac- 
terized by a sharp drop at moderate surface brightness 
values, followed by a flattening in the high surface bright- 
ness tail. This bi-modality reflects the transition from 
a regime in which the flux is dominated by the contin- 
uum contribution (the low surface brightness end) to a 
regime dominated by the emission of bright lines (the 
high surface brightness end). The transition occurs at a 
surface brightness of /pc ~ 20 ph s _1 cm~ 2 sr^ 1 keV^ 1 , 
which therefore we take as the reference value for the 
Galactic foreground continuum signal in this paper. The 
other curves of Fig. [3] refer to different energy resolu- 
tions. As the energy resolution degrades, the fraction 
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of high Galactic surface brightness becomes larger. The 
influence is considered in § 16.21 

Having estimated the two noise sources (/cxb and /fg 
to 30 and 20 ph s _1 enr 2 sr _1 keV -1 , respectively) and 
ignoring the negligible cont ribution from the instrumen- 
tal noise (|Kellev et al. 2007), one can compute the signif- 
icance of a line detection from the following expression: 



g ' inC ^[/nne + (/cXB + / FG )Ai ; ;] V ""^ ff ' W 

where /n n e is the source line flux in units of ph s _1 cm -2 
sr^ 1 , AE is the energy resolution of the instrument, A e g 
its effective area, AD its angular size (or the pixel size 
in the case of this work), and t oxp is the exposure time. 
Eq. ([]} assumes that both signal and noise can be treated 
as Poisson variables. 

In Table Q] we list for the two WHIM model explored 
(column 1) and two different exposure times considered 
in this work (column 2), the minimum line surface bright- 
ness required for a 5a line detection calculated with 
Eq. (JTJ) in a hypothetical observation performed with 
CRIS (col. 3). The corresponding number of detections 
per resolution element and unit redshift can be obtained 
from the dN/dz statistics plotted in Fig. [2j They are 
listed in columns 4 (O VII), 5 (O VIII), and 6 (simulta- 
neous O VII and O VIII detections). Column 7 shows the 
number of O VII+O VIII WHIM detections expected in 
the 1 degree 2 FOV due to the gas at z < 0.5. Each entry 
in columns 4-6 shows two values. The higher one repre- 
sents the expected number of detected lines contributed 
by all gas, not just in the WHIM phase. The smaller 
one, in parenthesis, refers to lines that can be attributed 
to the WHIM. All detection estimates assume an angular 
resolution of 2.6' x 2.6' (corresponding to the sum of 2x2 
pixels in the simulated dataset to the CRIS goal resolu- 
tion) and 1 eV energy resolution. These estimates do not 
account for the presence of spurious emission lines that 
contaminates ^20% of the energy band analyzed. Their 
impact on the ideal detection estimates listed in Table [T] 
is estimated in § [6j 

^From Tabled] we see that the expected number of de- 
tections is quite large even in the least favourable case 
(model Bl, i cxp =100 ks, simultaneous O VII+O VIII de- 
tection). However, only 1 Ms long observations guaran- 
tee to observe WHIM lines in the most conservative sce- 
nario of model Bl. With the more realistic model B2 one 
expects to be able to detect WHIM lines also with shorter 
(100 ks) exposure times. We note that WHIM detection 
estimates in column 7 are to be regarded as conserva- 
tive since they only consider gas out to z = 0.5 whereas 
hydrodynamical simulations show that the WHIM mass 
fraction is still quite large at z = 1.0. Overall, these 
results indicate that next-generation instruments will al- 
low one to unambiguously detect the WHIM in emission 
and, as we will see in the next sections, to characterize 
its thermal state and trace its spatial distribution. 

4. PHYSICAL PROPERTIES OF THE WHIM IN EMISSION 

In this section we investigate the physical properties 
(density p gas , temperature T, and mass M) of the WHIM 
that, according to our results, can potentially be detected 
by deep resolved spectroscopy in next-generation X-ray 
missions. The relevant quantities obtained from the anal- 



ysis of the mock X-ray spectra will be compared with 
the corresponding "true (input)" physical properties of 
the WHIM measured directly from the simulation. To 
perform realistic comparisons the "true" quantities were 
measured by averaging over all gas particles that con- 
tribute to the spectral signal, i.e. particles within vol- 
ume elements. The data used for comparison are taken 
from the FOV of 5.5° x 5.5° and the redshift range of 
0.202 < z < 0.274. The angular size of each volume ele- 
ment is 2.6' x 2.6' (the angular resolution) in § 14.11 and 
1.3' x 1.3' in § 14.21 The depth of a volume element is 3 
Mpc (comparable to the energy resolution). 

4.1. Density and temperature of detectable WHIM 

The thermal state of the gas can be appreciated from 
the phase-space diagram shown in top-left panel of Fig.UJ 
Colour contours are drawn in correspondence of the same 
gas mass fraction, indicated in the color -scale. This plot 
is ana logous to the one shown in Fig. 6 of lBranchini et ahl 
(pOOl with some important differences. First of all here 
we only consider particles with T > 10 5 K that are 
potentially relevant for the line emission. Second, this 
plot refers to redshift slice 0.202 < z < 0.274, consid- 
er ably narrower than < z < 0.2, the one considered 
bv lBranchini et ail (|2009ft . Finally, and most important, 
the plot refers to gas density and temperature averaged 
in the 2.6' x 2.6' x 3 Mpc volumes and not to the density 
and temperature of each gas particle in the hydrody- 
namical simulation, which cannot be observed directly 
The simulated data considered in Fig. |4] refers to the 
full FOV of 5.5° x 5.5° but rather thin redshift slice to 
avoid large variation in the volume of the resolution el- 
ements. The smoothing effect induced by the volume 
averaging is very evident in the poor sampling of the 
high-density and high-temperature regi ons, when com- 
pared to Fig. 6 of lBranchini et afl <|2009f ) . The effect can 
be further appreciated by looking at the black curve in 
the bottom panel of Fig. 0] which shows the probability 
distribution function of the gas density mass-weighted 
within each resolution element. The gas density on the 
X-axis is expressed in units of mean density. The dis- 
tribution function, which was obtained by marking the 
phase-space plot over the gas temperature, looks remark- 
ably symmetric. This symmetry is a result of the av- 
eraging process that deletes the high density tail that 
character ize the density di stribution function of the gas 
particles (jDave et al.1120011 ). 

The top right panel of Fig. 0] is analogous to the top 
left one but considers only gas elements for which the 
combined O VII and O VIII line surface brightness is 
above the 5 a detection threshold of 0.07 ph s _1 cm~ 2 
sr _1 , line systems that could be detected with 1 Ms 
observation with CRIS. This plot clearly indicates that 
next-generation instruments will preferentially probe the 
hottest part of the WHIM with T > 10 6 K with typ- 
ical density S a ~ 1 00, hence confirming the results of 
iBertone et alT(|2OT0() . 

This result is better illustrated by the red dotted curve 
in the bottom panel of Fig. |4] which is analogous to the 
black curve but only refer to the gas that can be poten- 
tially detected in a 1 Ms observation. It clearly shows 
that detectable emission lines preferentially probe high 
density regions. The integral of the the red dotted curve 
shows that 30% of the total gas mass (the integral of 
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Fig. 4. — Top left: Contours of constant gas mass fraction in phase-space, i.e., the fraction in unit logp and logT intervals. The gas 
density is normalized to its cosmic mean (X— axis) and its temperature is in unit of K (V-axis). The plots considers only gas particles in 
the redshift slice 0.202 < z < 0.274 and with temperature T > 10 5 K. Color coded contours are drawn in correspondence of different values 
of the gas mass fraction, indicated in the color scale. Top right: same as the top-left panel, but referring to gas elements characterized by 

VII + O VIII line systems strong enough to be detected with 1 Ms observation with CRIS. Bottom panel: the probability distribution 
function of the gas density obtained by merging the gas mass fraction in the p-T plane over its temperature. Black curve: all gas elements. 
Red dotted curve: gas elements in with both O VII and O VIII above 5 <r detection threshold of 0.07 ph s _1 cm -2 sr _1 , corresponding to 

1 Ms exposure with CRIS. Blue dashed curve: same as the red dotted curve but referring to a 100 ks exposure with detection threshold of 
0.48 ph s _1 cm -2 sr _1 . 



the black curve) could be probed via emission line spec- 
troscopy. The remaining 70 % of the gas that will go 
undetected typically resides in regions with density be- 
low 40 times the cosmic mean and constitutes the bulk 
of the WHIM. Reducing the exposure time to 100 ks 
further decreases the fraction of detectable gas to ~ 20% 
and shifts its mean density to larger values (blue dashed). 

To further investigate the issue of the WHIM de- 
tectability, we plot in Fig. [5] the contours of gas mass 
fraction as a function of gas density (X-axis) and line 
surface brightness (K-axis). The two plots refer to the 
O VII triplet (left panel) and O VIII line (right panel), re- 
spectively, and consider all gas elements with T > 10 5 K 
as in the top-right panel of Fig. |4j The first thing to 
notice is the very large range (about 10 order of mag- 



nitudes) of line surface brightness. Deep observations 
with a Xenia like mission will be capable of probing the 
gas above a detection limit of ~ 0.1 ph s _1 cm~ 2 sr _1 . 
The expected line surface brightness spans ^2 order of 
magnitude when S g > 100. Below this value, and for 
a line surface brightness below the detection threshold, 
the iso-probability contours run almost parallel to the 
K-axis, indicating the difficulty of detecting gas in lower 
density environment: huge improvements in the instru- 
mental sensitivity would be required to probe the bulk 
of the WHIM in emission. 

4.2. Estimating the Gas Temperature 

The quantitative analysis of the mock spectra allows 
us to investigate the thermal properties of the emitting 
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Fig. 5. — Contours of constant gas mass fraction in the line surface brightnesses-gas density plane. The plots considers only gas particles 
in the redshift slice 0.202 < z < 0.274 and with temperature T > 10 5 K. The surface brightness (Y-axis) is in units of ph s _1 cm" 



and the gas density (X-axis) is normalized to its cosmic mean. Left panel: Surface brightness of the O VII triplet. Right panel 
recombination line. The detection limit with a Xenia 1 Ms observation is 0.07 ph s _1 cm~ 2 sr -1 , or Y = —1.15. 
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Temperature measured from OVIII/OVII ratio (10 6 K) 

Fig. 6. — Relation between the gas temperature measured from 
the O VIII-O VII line intensity ratio in the mock WHIM spectra 
(X— axis) and different estimators for the underlying temperature 
of the gas (Y-axis). The Y— axis of the black, red, and blue crosses 
refer to O VH-cmission weighted temperature, O Vlll-cmission 
weighted temperature, and density weighted temperature, respec- 
tively. The data are taken from the FOV of 5.5° X 5.5° and the 
redshift range of 0.202 < z < 0.274. The angular size of each vol- 
ume element 1.3' X 1.3'. Horizontal error bars correspond to the 
temperature bin size, while the vertical bars represent the scatter 
among resolution elements in the temperature bin. The perfect cor- 
relation (i.e., Y = X) is represented by the central solid line. The 
two external dashed lines represents Y = l.ZX and Y = X/1.3, 
respectively. 

gas, namely to estimate the gas temperature. For those 
spectra in which both the O VII and O VIII lines have 
been detected one can use the line ratio to estimate the 
gas temperature. Alternatively one could estimate the 
gas temperature by comparing the relative intensity of 
the resonance, forbidden, and intercombination lines of 
the O VII triplet. These, however, can be detected only 
in the ~ 20 % of the cases. Therefore in this section 
we assess how well one can measure the gas temperature 
from the O VII-0 VIII line intensity ratio when the CIE 
hypothesis holds true. For this purpose we have com- 
pared the temperature estimated from the mock spectra 



to the "true" temperatures of the gas responsible for line 
emission. The latter quantity, however, is not uniquely 
defined for two reasons. First, there is no one-to-one cor- 
relation between the measured quantity, which is derived 
from the line intensity, and the temperature of the gas 
particles. Second, the line emitting regions are sampled 
in finite volume elements because of the angular /energy 
resolution of the spectrometer. To tackle this problem 
we will use different estimators of the "true" temperature 
and compare the results with the temperature estimated 
from the O VII-0 VIII line intensity ratio. 

In Fig. [6] we have plotted the correlation between 
the temperature measured from the mock spectra (X— 
axis) and the one obtained from three different gas 
temperature estimators (Y-axis). Here we consider 
the case of model B2 only. In the plot the dif- 
ferent colors indicate the three different estimators: 
O VH-line-emission weighted temperature (black crosses) 
O VHI-line-emission weighted temperature (red crosses) 
mass-density weighted temperature (blue crosses). The 
"true" gas temperature is estimated in those resolution 
elements where the surface brightness of both O VII and 
O VIII is above 0.07 ph s _1 cm -2 sr -1 . The vertical size 
of each cross represents the scatter about the mean. The 
horizontal bar represents the size of the corresponding 
temperature bin. 

The plot shows that a tight correlation between the 
measured and the "true" temperature is found with all 
estimators used. The temperature estimated from the 
O VIII-0 VII line intensity ratio seems to best trace the 
O VH-emission weighted temperature of the gas, which is 
expected since temperature estimated from spectral lines 
preferentially weight those regions in which line emission 
is stronger. When the O VHI-emission weighted temper- 
ature is used the correlation is still tight but the slope is 
larger than unity. This is due to the fact that this estima- 
tor preferentially weights regions in which the O VIII is 
strong, i.e. regions in which the temperature is typically 
higher than those from which the O VII photons come 
from. With the mass-density weighted temperature the 
relation is unbiased for temperatures below T ~ 3x 10 6 K 
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and then the correlation becomes steeper. This behavior 
can be understood as follows. The temperature mea- 
sured from the spectra is sensitive to T =(2 -3)xl0 6 K. 
where O VII and O VIII emissivity has a peak and it 
is reliably traced by the mass weighted temperature. At 
higher temperatures, a large fraction of the gas is too hot 
to contribute to the O VII emission and therefore it is not 
traced by the O VII line. As a result, the measured tem- 
perature underestimates the mass-weighted temperature 
of the gas. 

Our assumption of CIE is well justified for gas in high 
density and temperature environments (8 g > 100 and 
T > 2 x 10 6 K). Photoionization becomes non- negligible 
in regions characterized by lower density and tempera- 
ture (i.e. when 5 g < 100 and T < 2 x 10 6 K). The pres- 
ence of the photoionizing background would change the 
ionization state of the gas. In par ticular, this backg round 
would ionize O VII into O VIII (jChen et alJl2003f >. and 
thus enhancing the intensity of the O VIII emission line 
relative to that of O VII. As a consequence, assuming 
the O VII/O VIII ratio of CIE overestimates the temper- 
ature of the gas. However, these effects will affect the gas 
properties in regions of moderate overdensity and tem- 
perature which could hardly be observed through emis- 
sion line spectroscopy with current and next-generation 
detectors. 

These results show that one can use the measured 
O VII and O VIII lines to estimate the temperature of 
the gas with a reasonable precision (10-20% uncertainty, 
depending on the gas temperature). The measured quan- 
tity probes the temperature of the gas responsible for 
O VII emission and underestimates that of the O VIII 
emitting regions. The measured quantity also traces the 
temperature of the bulk of the gas for T =(1-3) xlO 6 K 
and underestimates it in hotter regions. 

Finally, we note that we have assumed a detector with 
angular resolution 1.3' x 1.3' and that we do not con- 
sider the line broadening due to the detector response. 
The rationale behind this choice of the angular resolution 
is to reduce possible systematics in the estimate of the 
"true" temperature. In fact, using an angular resolution 
of 2.6' x 2.6' as in the rest of the paper would not change 
appreciably the value of the measured temperature but 
would artificially reduce the estimated "true" tempera- 
ture of the gas because of the increase in the filtering 
volume associated to the resolution element. A good an- 
gular resolution is thus important for an accurate deter- 
mination of the physical state of the WHIM. Including 
the effect of the line broadening induced by the energy 
resolution of the detector would reduce the number of 
detectable line systems and would increase the statisti- 
cal errors. However, we do not expect that these errors 
would seriously affect our temperature estimates since 
temperature errors are driven by the presence of multi- 
phase gas within the resolution elements and not by the 
limited line statistics. 

5. PROBING THE SPATIAL DISTRIBUTION OF THE WHIM 

Spatially resolved spectroscopy provides the unique 
possibility of tracing the spatial distribution of the 
WHIM. In this section we use our simulated 2D spec- 
tra to investigate the 3D distribution of the WHIM re- 
sponsible for detectable line emission. Here we focus our 
attention on the WHIM model B2 and assume a 1 Ms 



observation with a CRIS-like spectrometer with angular 
resolution of 2.6' x 2.6', energy resolution of 1 eV and 
effective collecting area of 1000 cm 2 . 

We only consider systems in which both O VII and 
O VIII lines can be detected at > 5tr significance and 
ignore the possibility of using additional metal lines as- 
sociated to the source. The need of using both lines is 
dictated by the need of measuring the redshift of the asso- 
ciated line emitting region. This is guaranteed by select- 
ing line pairs with energy ratio consistent with that of the 
O VII and O VIII lines. Using additional metal lines can, 
in principle, improve the identification of line emission 
systems and the determination of their redshift. How- 
ever, the benefit in using these extra lines is in practice 
limited by the associated increase in chance contamina- 
tions. For example, a number of emission lines from Fe L 
complex can be seen in the mock spectra. However, these 
lines typically probe the hotter part of the WHIM, typi- 
cally associated with virialized structure also responsible 
for strong continuum emission which reduces the signif- 
icance of line detections. The Ne IX is a better signpost 
for the WHIM since it probes gas with T ~ 10 6 6 K. 
However its emission line is too close to the Fe L lines to 
eliminate the risk of spurious detection. Finally, the en- 
ergies of the redshifted N and C lines are typically found 
at the edges of the instrumental range, where a forest 
of emission lines due to L and M transition of heavier 
metals exists. For all these reasons we search the spectra 
for emission systems characterized by the O VII Ka reso- 
nance and O VIII Lya pair and do not look for additional 
ion lines. 

Fig. [7] and [H] compare the 3D distribution of the 
selected O VII + O VIII line systems (lower panels) 
with the spatial distribution of the gas in the simula- 
tion (upper panels). Both figures show the same FOV of 
5.5° x 5.5°, corresponding to a 6 x 6 mosaic with CRIS, 
but refer to two different redshift slice. Fig. [7] refers to 
the redshift interval z = [0.2117, 0.2317] corresponding to 
a 54/i _1 Mpc slice at a comoving distance of 630 /i -1 Mpc 
in the cosmological model of the hydrodynamical simu- 
lation. Fig. [5] shows a much closer slice of 57/i _1 Mpc at 
a comoving distance of 265 /i _1 Mpc, corresponding to 
z = [0.0805,0.1004]. 

In the gas map the blue regions correspond to iso- 
density surfaces drawn at p gas / (Pgas) = 75. The green 
surfaces corresponds to regions characterized by smaller 
overdensity /9 gas / (Pgas) = 10. At these overdensity lev- 
els the spatial distribution of the gas in the simula- 
tion is characterized by the well-known network of fil- 
aments punctuated by virialized structures that, in this 
iso-contour maps, appear as quasi-spherical "blobs" . The 
filamentary structure is more clearly seen in the more dis- 
tant slice thanks to the large volume sampled. The gas 
distribution in the nearby slice is clearly more dominated 
by a few, prominent structures, as expected. 

The red regions in the bottom panels identify the posi- 
tions of the selected O VII + O VIII line systems. Clearly 
these line systems qualitatively trace, in a sparse fash- 
ion, the large scale structures of the gas with overden- 
sity Pgas/ (Pgas) = 75. However the tracing is not un- 
biased since the line systems tend to undersample the 
more massive virial structures associated to the largest 
blob. This bias reflects the unfavorable ion balance for 
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Fig. 7. — Slice of gas density contour (top) and O emission 
map (bottom) at z = 0.2117 — 0.2317, corresponding to comoving 
distance of 603-657/1- 1 Mpc. The field of view is 5.5° X 5.5°. The 
blue region of the top panel corresponds to p ga s/ (Pgas) = 75, while 
the green region corresponds to Pgas/ (Pgas) = 10. The red regions 
(blobs) show where both O VII and O VIII lines are identified 
with > 5cr significance, with 1 Ms exposure by Xenia CRIS-like 
instrument. 

O VII and O VIII within hot virialized regions. The fact 
that O VII+O VIII line systems are poor tracers of gas 
at lower overdensity (/c ga s/ (Pgas) = 10) confirms the fact 
that emission lines preferentially samples the fraction of 
the WHIM that resides in high density regions (Fig.@|. 

This qualitative analysis confirms that next-generation 
X-ray spectrometers could be successfully used to trace 
the spatial distribution of the WHIM through its line 
emission. Thanks to the large number of the de- 
tected emission line systems, more quantitative analy- 
sis can be performed through statistical estimators. In 




Fig. 8.— Same as Figured but slice at z = 0.0805 - 0.1004, 
corresponding to comoving distance of 237-294/i — 1 Mpc. The field 
of view is 5.5° x 5.5°. 

lUrsino et all (|2010aD . which uses the same mock spectra, 
it is shown that the two-point spatial correlation function 
of the line-emitters could be measured with a good ac- 
curacy. Its slope and correlation length, which should be 
similar to those of the galaxies, could be estimated at a 
few % accuracy level. 

The tracing of underlying gas by the detected line sys- 
tems, which in the plots looks quite patchy, is possibly 
improved by adopting some adaptive smoothing tech- 
nique like those based on the Wavelet transform anal- 
ysis. The improvement would be particularly noticeable 
for nearby structures since, in this cases, the angular size 
of the associated gas cloud is larger than the angular res- 
olution of the instrument. 
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6. ROBUSTNESS TESTS AGAINST NON-IDEAL EFFECTS 

In this section we assess the impact of foreground emis- 
sion, energy resolution and contamination from spurious 
metal lines on the detectability of the extragalactic O VII 
and O VIII lines and on our capability of tracing the spa- 
tial distribution of the WHIM. 

6.1. Impact of contamination from spurious lines and 
high continuum 

We address the problem of identifying O VII and 
O VIII lines with no a priori knowledge of the nature of 
the detected lines and of their redshift. For this purpose 
the line-detection procedure have been modified from 
that depicted in § 13.21 The rationale behind this choice 
is to avoid contamination from lines associated to differ- 
ent line emitting regions along the line of sight and to 
account for the fact that the line detection efficiency is 
not constant across the FOV since the level of the contin- 
uum emission depends on the direction. We have decided 
to ignore additional metal lines associated to the WHIM 
since their search would increase the chance of false iden- 
tifications. To identify the putative WHIM emission sys- 
tems we do not perform spectral fitting assuming some 
particular model but we search for the O VII Ka and 

VIII Lya line pairs with a strategy more conserva- 
tive than that adopted in § [3] and designed to account 
for the variable spectral continuum and for the presence 
of spurious spectral features. The procedure is as fol- 
lows. (1) We deal with the fact that the continuum level 
may vary along the spectrum and also from spectra to 
spectra, across the FOV. To account for this effect we 
modify the line detection criterion: instead of setting a 
fixed line detection threshold, corresponding to a mini- 
mum line surface brightness, we define the detection sig- 
nificance S by comparing the number of photon counts in 

1 eV bin ( Counts) with the number of continuum counts 
[Continuum) obtained by averaging over the neighbor- 
ing 50 cV that do not contain line features above 2 a or 
known galactic lines brighter than > 0.05 ph s _1 cm" 2 
sr" 1 eV" 1 . The significance is then defined as: 

Counts — Continuum , . 

S/N = (2) 

1.0 + v Continuum + 0.75 

whe re the denomin ator is an approximate la Poisson er- 
ror (lGehrelslll986lh Note that Eq. © is a similar quan- 
tity as Eq. ([1} but is a quantity calculated with observed 
quantities. (2) We consider only lines in the energy range 
[375,653] eV detected at 5a significance (i.e., S/N > 5). 
With this criterion we identify 7.8 x 10 3 lines deg -2 corre- 
sponding to transitions of highly ionized ions of N, O, Ne, 
and Fe. (3) We minimize the chance of contamination by 
spurious lines by enforcing a number of constraints, (i) 
we reject all lines with flux of > 0.05 ph s" 1 cm" 2 sr -1 
eV" 1 that could be attributed to Galaxy emission, i.e., 
all known Galactic lines; (ii) to avoid oversampling the 
same structures we discard all lines with energy within 
15 eV from that of a stronger emission line; (Hi) we search 
for line pairs from the same angular element, and pick 
up the pairs only if the energy interval of the two lines 
is the same as that of rcdshiftcd O VII Ka resonance 
and O VIII Lya lines. The emission lines without an 
associated O VII or O VIII line are discarded. 



Using the procedure, we only consider those systems 
that both the O VII Ka and the O VIII Lya lines are 
identified with significance larger than 5 a. As a result we 
obtain 446 deg - line emission systems identified. This 
number is ~ 30% less than those identified in § O The 
discrepancy comes from the combined effect of line con- 
tamination and high continuum level. We confirm that, 
with this procedure, the contamination of chance coinci- 
dence of either O VII or O VIII lines with heavier metal 
lines is less than 10%. Our spectra do not include pho- 
tons emitted by gas at z > 0.5. O VII lines emitted 
beyond this redshift would be redshifted out of the en- 
ergy range probed by the detector, so the only impact 
of high redshift gas would be to increase the probability 
of contamination by higher energy lines such as O VIII, 
Ne IX, O X and Fe L complex. However, we do not ex- 
pect a severe additional contamination for two reasons. 
First, the dimming of the emission line intensity, which 
scales as (1 + z)" 3 , is not exactly counterbalanced by the 
increase in the probed volume since, the latter decreases 
more slowly than cx (1 + z) 3 at high redshifts Second, at 
higher redshifts the mass fraction of the gas in the Ly- 
a forest increases while the abundance of the gas con- 
tributing to O VII and O VIII lines emission decreases. 
To test the validity of these hypotheses and quantify line 
contamination by high-redshift gas we have created and 
analyzed mock spectra that accounts for emission by the 
gas out to z < 1.2 and within a FOV of 2.6° x 2.6° 
and have compared the results with the reference case of 
z < 0.5. Including the emission from gas beyond z = 0.5 
decreases the number of detectable O VII-0 VIII pairs 
by only 10-15%.. 

6.2. Impact of energy resolution 

A degradation in the energy resolution of the instru- 
ments affects the possibility of detecting the WHIM lines 
in two ways. First, it decreases the signal-to-noise ratio 
in Eq. fl}. As a consequence, one needs to increase the 
surface brightness threshold for 5a detections, hence re- 
ducing the total number of detected lines. Second, with a 
large AE the chance of contamination from neighboring 
lines increases. Here we let the energy resolution vary 
in the range [1,7] eV which encompasses the values al- 
ready ac hieved with current instrument (AE = 7 eV with 
Suzaku; iKellev et al.l 12007ft and the energy resolutions 
expected by currently planned X-ray missions equipped 
with a microcalorimeter spectrograph. 

A visual impression of the cumulative impact of these 
two effects is given by Fig. |9l which shows a portion of 
the mock emission spectrum displayed in Fig. [T]but ob- 
served with different energy resolutions and in log scale. 
With 1 eV resolution (FWHM), O VII triplets and O VIII 
resonant line at z = 0.033 (0.55 keV and 0.63 keV, re- 
spectively) are clearly detected, with a marginal detec- 
tion of a emission line at 0.47 keV. With AE = 3 eV, 
the 0.47 keV line is below detection threshold, and with 
AE = 5 eV the O VII triplet becomes hard to observe 
because of the local (z = 0) O VII emission. When AE is 
7 eV, the O VII triplet drops below detection threshold. 

A quantitative assessment of the contamination by 
spurious lines is provided by Fig. [3] The fraction of 
energy bins with surface brightness above that of the 
Galactic foreground emission (/fg = 20 ph s" 1 cm" 2 
sr" 1 keV" 1 ) steadily increases with AE: from 20% for 
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TABLE 2 

Number of emission line detections with different detector energy resolution 



AE / linc <2iVoviI+OVlll/<fe Fraction of Novii+OVIII per deg 2 

(eV) (ph s^ 1 cm - 2 sr^ 1 ) ■ • • high /pq (%) in the simulation 

— a 0.07 2.4 14 446 

1 0.07 2.4 20 396 

3 0.11 2.1 30 351 

5 0.13 2.0 36 274 

7 0.15 1.9 41 163 



Note. Column 1: Detector energy resolution. Column 2: Minimum line 



surface brightness required for a 5a detection (ph s _1 cm -2 sr" 1 ). Column 3: Ex- 
pected number of simultaneous O VIII and O VIII detections per resolution element 
and unit redshift. The foreground emission is assumed to have a constant surface 
brightness of 20 ph s _1 cm -2 sr _1 keV -1 . Column 4: Fraction of the energy range 
where foreground emission surface brightness exceeds 20 ph s _1 cm~ 2 sr _1 keV -1 . 
Column 5: Number of simultaneous O VIII and O VIII detections per square degree 
in the mock spectra. All estimates assume an angular resolution of 2.6' x 2.6' and 
model B2. 

a Not convolved with detector energy resolution. The bin size is 1 eV. 

< 20%. On the other hand, contamination by galactic 
foreground has a more serious impact since the num- 
ber of expected detection is directly proportional to that 
of the uncontaminated energy bins. Our results show 
that decreasing AE from 1 to 7 eV reduces the number 
of expected detection by a factor of 2. It should also 
be noticed that this effect does not randomly affect all 
emission line systems. The systems that are more se- 
riously affected are those with redshifted emission lines 
that coincide with those of the Galactic foreground. This 
induces a non-trivial selection effect in the redshift distri- 
bution of the WHIM emitters that need to be corrected 
for to trace their 3D distribution. 

To quantify the impact of the energy resolution on the 
expected number of WHIM line detections we have re- 
peated the analysis performed in § 16.11 after convolving 
the mock spectra with the appropriate response matrix. 
In the analysis of the response-convolved spectra, we 
again do not fit the spectra with a particular model, but 
searched for the O VII and O VIII line pairs. The only 
difference is that in the new analysis photon counts in 
neighboring energy bins are summed up before search- 
ing for emission lines. This is done to collect photons 
that are spread into several energy bins (the bin width 
is 1 eV) due to the energy resolution of the instrument. 
The results of this procedure obviously depend on the 
number of neighboring bins to be summed upon. We 
set this free parameter by maximizing the number of de- 
tectable emission lines. It turns out that the best results 
are obtained when summing over {{AE /I eV) + 1) bins 
for AE < 5 eV and over {{AE/1 eV) -3) for AE = 7 eV. 
The change of the trend for AE = 7 eV reflects the fact 
that at low energy resolution foreground contamination 
becomes the limiting factor. The number of expected de- 
tections is listed in column 5 of Table [5J It is clear that 
an energy resolution AE < 3 eV is desirable, although 
even with AE = 7 eV one still expects 160 O VII+O VIII 
line detections per square degree. The reason why the 
number of expected detection of the 'no convolution' case 
is different from the analogous one listed in Table [1] is be- 
cause in the line detection strategy used in this case, and 
described in § 16.11 we have ignored the energy range in 
which the surface brightness of the foreground emission 




Energy (keV) 



Fig. 9. — Mock spectrum in a 2.6' X 2.6' area from a 1 Ms ex- 
posure with Xenia CMS. The black spectra are the sum of Galac- 
tic foreground and unresolved extragalactic background emission, 
while the red is extragalactic diffuse gas. The four panels of the 
each plot show the spectra convolved with the energy resolution 
(FWHM) of 1 eV, 3 eV, 5 eV, and 7 eV, respectively. 

AE = 1 eV to 41% for AE = 7 eV. In Tabled] we quan- 
tify how energy resolution affects the number of expected 
detections, splitting its effect in two: the increase of the 
surface brightness detection threshold (column 3) and 
the increase in the number of bins contaminated by the 
Galactic Foreground (column 4). All estimates assume 
the same setup as in Tabic Q] apart from energy resolu- 
tion: WHIM model B2, 1 Ms observation, and 2.6' x 2.6' 
angular resolution. The first effect is rather minor and 
reduced the expected number of linc detection only by 
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is above /fg- These ranges could be also be searched 
for detectable lines by locally increasing the detection 
threshold. However the expected increase in the number 
of detections is small since the intensity of the foreground 
emission lines is typically an order of magnitude larger 
than that of the WHIM lines. 

7. DISCRIMINATING AMONG WHIM MODELS 
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Fig. 10. — The expected dN/dz curve with a Poisson error with 
1 Ms exposure, 2.6' X 2.6' angular size, and 1 deg 2 FOV with CRIS. 
The black curve is the model curve (model B2). Color-filled areas 
arc la Poisson error regions: cyan and blue are without considering 
detector energy resolution and with convolution with 7 eV energy 
resolution, respectively. Red dashed line is from model Bl (Poisson 
error is not shown). 

As discussed in § [2l current limitations in our under- 
standing of the stellar feedback, chemical evolution, and 
metal diffusion process are reflected in the large uncer- 
tainties of the WHIM models. The most recent hydro- 
dynamical simulations of the intcrgalactic gas at low red- 
shifts have been used to construct different WHIM mod- 
els by assuming different, physically plausible, recipes 
for the above processes. Comparisons between model 
predictions and observations can therefore constrain the 
cosmic history of chemical evolution and metal diffusion 
along with the underlying stellar feedback processes by 
discriminating among these WHIM models. UV obser- 
vations are leading the way, since they currently pro- 
vide the most stringent constraints on the WHIM mod- 
els through O VI line statistics. X-ray observations can 
add additional constraints, through line statistics, direct 
assessment on the physical state and spatial distribution 
of the emitting gas. 

In this section we assess to what extent future exper- 
iments will enable us to discriminate among different 
WHIM models. Here we focus on the dN/dz statistics of 

VII and O VIII lines. Fig. ITOl shows the expected num- 
ber of O VII+O VIII detection per unit redshift in an an- 
gular resolution element of 2.6' x 2.6' and for 1 Ms obser- 
vation with CRIS for model B2 (black, continuous curve) 
and Bl (red, dashed). The cyan-colored band around the 
black lines show the la Poisson uncertainties for model 
B2 from the finite number of detections expected in a 

1 deg 2 FOV. The expected errors are much smaller than 
the difference between model Bl and B2, suggesting that 



observations will efficiently discriminate among different 
WHIM models. The capability of constraining the phys- 
ical processes that regulate the WHIM is further illus- 
trated by pointing out that difference between model Bl 
and model B2 is of the sam e order as that betw een the 
WHIM models proposed bv iCen fc Fanjl (|2006[ L which 
adopt different prescription for the ionization balance 
and for galactic superwinds. The larger, blue uncertainty 
strips that extend beyond the surface brightness thresh- 
old of 0.15 ph s _1 enr 2 sr -1 refers to the case of AE = 7 
eV. It shows that even with a rather modest energy reso- 
lution the discriminatory power of the dN / dz test is still 
considerably large. The significance of discriminating the 
two models is 3-5cr with a 1 deg 2 observation, depending 
on the energy resolution. 

It should be noted that the cosmic variance introduces 
an additional error in determining the dN/dz curve. 
We have estimated the contribution of cosmic variance 
by splitting the simulated 5.5° x 5.5° field into several 
smaller regions. The variance in detectable emission line 
systems is 20%, 15%, and 10%, for the FOV of 0.25 deg 2 , 
1.0 deg 2 , and 4.0 deg 2 , respectively. Hence, a large- field 
mapping with > 4.0 deg 2 is desirable to make the uncer- 
tainty due to cosmic variance smaller than the statistical 
error (^10% for 7 eV energy resolution). This can be 
easily achieved with a dedicated mission such as EDGE 
and Xenia thanks to the large FOV. 

Finally, we point out that emission-line studies provide 
better statistics than absorption studies because of the 
comparatively larger number of line detections. In fact, 
the two approaches are very complementary. Emission 
studies preferentially probe regions of enhanced density 
where the impact of stellar feedback is larger and chemi- 
cal enrichment is presumably stronger. On the contrary, 
absorption studies preferentially probe regions of mod- 
erate overdensity where the bulk of the WHIM resides 
and thus are best indicated for closing the baryon bud- 
get at z = 0. Space missions like EDGE, Xenia, and 
ORIGIN are designed to perform both emission and ab- 
sorption analyses thanks to their capability of fast re- 
pointing that will enable us to use Gamma-ray bursts as 
distant beacon and search for the chara cteristic absorp- 
tion features of the intervening WHIM (jBranchini et al.l 
[2001 . 

Constraining abundance and physical properties of the 
WHIM from observations it is not straightforward. The 
first reason is that in order to trace the thermal his- 
tory of the WHIM it requires an accurate and self- 
consistent treatment of metal line cooling which, instead, 
is ignored in the present work. The second reason is 
the approximate treatment of a number of processes 
and quantities (galactic winds, AGN fe edback, stellar 
IMF and sub-resolution turb ulence; see iBertone et al.1 
l20Ul ?; ?: lSchave et 10120101 ) that may affect the ther- 
mal state of the WHIM in the current numerical ex- 
periments. Some of these processes have little impact 
on observ able quantities such as O VIII surface bright- 
ness (e.g.. IBertone et all l2010f ) . and hence it cannot be 
used to discriminate among competing WHIM models. 
On the other hand, emission line statistics and analysis 
will no doubt help to lift some of the model degener- 
acy, especially when coupled to observational constraints 
on the spatial distribution of the line-emitting gas, like 
those derived from the estimate of the angular and spa- 
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tial two-point correla tion function of the O-line emission 
(|Ursino et al.ll2010al) . 

8. DISCUSSION AND CONCLUSIONS 

In this work we have analyzed the mock X-ray spcc- 
tr a extracted from the the hydrodynamical simulation 
of iBorgani et al.l ((2004) to investigate the possibility of 
detecting the WHIM in emission with next-generation X- 
ray satellites, assess its thermal state and trace its spatial 
distribution. Our th eoretical predict i ons ar e based on 
the WHIM model of IBranchini et al.l (I2009I ). The the- 
oretical predictions presented in this work are mainly 
based on their model "B2" although we also consider 
the more conservative (and less realistic) WHIM model 
"Bl" to bracket theoretical uncertainties. Finally, the 
experimental setup considered here is the same one pro- 
posed for the Xenia mission. In particular we will assume 
that our mock spectra were "observed" with the imaging 
spectrometer CRIS. 

We showed that a CRIS-like spectrograph will enable 
us to detect a large number of O VII+O VIII line emit- 
ting systems, to estimate their temperature from the line 
ratio and to trace their spatial distribution. The line sys- 
tems are expected to trace the hotter and denser part of 
the WHIM, typically associated to the outskirts of large 
virialized structures. Our results suggest that the num- 
ber of expected detections is large enough to discrim- 
inate among different WHIM models and to constrain 
the physical processes that regulates its thermal state 
and metal content, even when one accounts for contami- 
nation effects and allows for a moderate (~7 eV) energy 
resolution. 

All the results presented in this work are based on a 
WHIM model that fulfill all observational constraints but 
is nevertheless poorly constrained by observations. From 
the theoretical side, WHIM models suffer from the lim- 
ited knowledge of the feedback mechanisms and metal 
diffusion processes that potentially affect the WHIM 
properties. In particular, the WHIM model adopted 
here, though effective in describing observations, lacks 
self-consistency since the metallicity of the gas is assigned 
in the post-processing phase rather than being computed 
during the run. To fix these problems one should mod- 
ify theoretical prescriptions for feedback process or metal 
diffusion and gas cooling. Let us asses qualitatively the 
impact of having ignored line cooling. Efficient metal 
cooling requires a cooling time t coo \ shorter than the Hub- 
ble time tn. With no metal cooling, i coo i < tn only when 
8 g > 1000. Cooling become more efficient in presence of 
metals and less efficient with the gas temperature, i.e. 
t CO o\ decreases with metallicity (e.g., Figure 4 of ?), and 
increases with T. On the other hand, we have shown 
that emission line spectroscopy can only detect gas with 
T > 10 6 K (Figure HJ). Taking all this into account and 
assuming a reference metallicity of ~ 0.1 Zq one finds 
that i cool < t H if S g > 300 for T = 10 6 K and if 5 g > 1000 
for T = 2 x 10 6 K. This simple arguments indicate that 
metal line cooling preferentially affects dense (6 g > 300) 
gas typically associated to galaxy groups rather than the 
WHIM. Indeed, inefficient cooling in the simulation is 
probably the cause of the excess X-ray emission from 
groups in the simulation. 

In fact, the impact of the metal cooling is likely to be 
more significant because it affects the whole thermal his- 



tory of the baryons (????). When metals are transferred 
from galaxies, either by galactic wind or AGN, both the 
density and metallicity are high, and hence metal cool- 
ing is very efficient. After the temperature is lowered by 
metal cooling, the gas is further cooled down by adiabatic 
expansion. As a result, the thermal state of baryons at 
z ~ is determined by their thermal history. The impact 
of these effects can only be assessed by running a numer- 
ical simulation in which metal cooling is accounted for 
self-consistently and, more importantly, in which the rel- 
evant physical process are included correctly. However, 
current numerical experiments whi ch include a self con- 
sistent treatment of metal cooling (|Bertone et al.ll2~010D 
seem to corroborate our qualitative argument since they 
show that metal cooling mostly affect regions with over- 
density typical of galaxy groups, rather than the WHIM. 
For this reason, we are quite confident that considering 
predictions from both model B2 and the very conserva- 
tive model Bl will allow us to bracket theoretical uncer- 
tainties, including incorrect treatment of metal cooling 
i n the numerical simul ation. 

iBertone et al.l (|201Q[ ) have shown that predictions on 
the detectability of the O VIII line are are quite ro- 
bust to the different model prescription with the excep- 
tion of the cooling function. They show that using a 
metal-dependent cooling instead of a metal-independent 
scheme (as in our model) decreases by up to an order of 
magnitude the SB of the O VIII line in regions populated 
by non- WHIM, but does not significantly alter the line 
emission properties in regions of lower overdensity where 
the WHIM typically resides. Reducing the artificially 
large brightness of large density regions would decrease 
the noise level and increase the probability of detecting 
the emission line systems associated to the WHIM. 

From these considerations, we conclude that our 
WHIM models are adequate to assess the detectability 
through O VII/O VIII line pairs of the WHIM in emis- 
sion and are confident that theoretical uncertainties are 
effectively bracketed, perhaps too generously, in the scat- 
ter among Bl and B2 model predictions. On the other 
hand, the gas distribution in the phase-space (Figure |4j 
and in the 3D space (Figure [7] and [8|), especially in re- 
gions of high density and moderate temperature, might 
not be accurate because of the approximate treatment of 
gas cooling in the simulations. 

The results presented in this paper assume collisional 
ionization equilibrium. This is a reasonable hypothesis 
since large part of O VII and O VIII lines that next- 
generation instruments are expected to detect are likely 
to be produced by warm-hot gas in high density envi- 
ronments where the CIE is a fairly good assumption 
(S g > 100 and T > 2 x 10 6 K). The validity of the 
CIE hypothesis can be directly checked by analyzing a 
number of spectral features. One is the relative com- 
parison of the intensity of the resonant, forbidden and 
intercombination lines of the O VII triplet, which also 
provide an estimate of the gas temperature, alternative 
to that obtained fro m the O VIII/O VII ratio (see, e.g., 
iPorquet et al.ll200ll ). The possibility of performing this 
test, however, is limited to the few cases in which the 
components of the triplets can be unambiguously de- 
tected. ^From the analysis of the mock spectra, only 
in 30 % of the cases can we detect both the resonant and 
the forbidden lines. The weaker intercombination lines 
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is only detected for ~ 10 %. If photo-ionization plays 
a role, O VIII line is also emitted from lower tempera- 
ture gas, in addition to O VII triplets. Therefore, the 
number of detected line systems would increase because 
lower density (S < 200) regions becomes also detectable. 
The ionization state of the gas in this low density re- 
gions would be affected and the temperature would be 
overestimated under the assumption of CIE. 

Our analysis of the spatial distribution of the line emis- 
sion systems presented in § [5] is rather qualitative. A 
more quantitative analysis, based on the mock spec- 
tra des cribed in this work, is presented in lUrsino et al.l 
(|2010af ) in which it is shown that the number density of 
detectable emission line systems is large enough to char- 
acterize the distribution by means of the spatial two- 
point correlation function. It is also shown that the 
two-point correlation function of the WHIM is signifi- 
cantly different from that of the non- WHIM and similar 
to that of galaxies, to which the line-emitting WHIM 
is associated. Its correlation properties do not seem to 
evolve significantly in the redshift range z = [0,0.5] and 
the dynamical state, induced from the velocity-driven 
anisotropics in the observed clustering, is consistent with 
that of a gas coherently falling toward larger structures. 

The main results of this work are summarized as fol- 
lows. While the analyses are based on the specification of 
CRIS onboard EDGE or Xenia, essentially the same con- 
clusions are applicable for DIOS or ORIGIN, which are 
equipped with a spectrometer of slightly different FOV, 
angular resolution, and effective area. 

• The WHIM can be detected in emission through 
O VII and O VIII lines. Unambiguous detection of 
the line systems requires the simultaneous detec- 
tion of both lines, each one with statistical signifi- 
cance above 5a. With a 1 Ms observation and in a 
FOV of 1 deg 2 we expect to detect - 640 O VII + 
O VIII line emitting systems. A significant number 
of detections (~ 300 per deg 2 ) is still expected with 
a much shorter (100 ks) observation. These results 
are in broad agreement with those of iCen fc FaneJ 
(l2006h . 

• A significant number of these line systems, how- 
ever, trace gas which is hotter or denser than the 
WHIM (5 > 1000). The number of line systems 
contributed by the WHIM is - 65 % (i.e., - 430 
deg -2 ) for a 1 Ms observation. The fraction drops 
to ~ 20 % when the exposure time is reduced to 
100 ks. In other words, deep observations are re- 
quired to detect a number of line systems large 
enough to investigate the statistical properties of 
the WHIM. Shorter observations would preferen- 
tially probe hotter gas in high density environments 
associated to clusters or groups. 

• The expected large number of detections is ro- 
bust. If one considers the very conservative model 
Bl, which we have implemented to assess theo- 
retical uncertainties, the number of line systems 
above 5a detection threshold is still large (~ 480 
deg~ 2 ), 40% of which contributed by the WHIM. 
The difference between model B2 and model Bl 
provides generous estimate of the theoretical uncer- 



tainties, significantly larg er than the scatter a mong 
the models predictions of ICen fc FansJ (|2006| ). 

• The number of detectable systems is affected by 
contamination from spurious lines and by the level 
of the continuum that may vary with energy and di- 
rection. It is confirmed that a more rigorous detec- 
tion procedure that accounts for both effects misses 
only ~ 30 % of the line systems that are originally 
selected. 

• The number of detectable emission lines decreases 
with the energy resolution. However, even with an 
energy resolution of 7 eV, already achieved by the 
microcalorimeter on board Suzaku, one expects to 
detect ~160 O VII+O VIII line systems per deg 2 
in 1 Ms observation. 

• Detectable line systems preferentially trace the 
denser regions of the WHIM with typical gas over- 
density of S g ~ 100, typically associated to virial- 
ized structures like galaxies, groups, and the out- 
skirts of galaxy clusters. Increasing the instrumen- 
tal sensitivity will not change these results signif- 
icantly since the bulk of the WHIM that resides 
in lower density environment is characterized by a 
very low line emissivity . Thes e results agree with 
those of iBertone et al.l (|2010D and illustrates the 
importance of studying the WHIM both in emis- 
sion (to assess its spatial distribution and thermal 
state) and in absorption (to probe low density re- 
gions and close the baryon budget). 

• In the region where CIE holds, the tempera- 
ture of the gas associated to the detectable emis- 
sion systems can be reliably estimated from the 
O VII/O VIII line ratio in the range of (1- 
5)xl0 6 K, with typical random errors of ~ 20 
%. The measured temperature trace the O VII- 
cmission weighted temperature of the gas, which, 
for T =(l-3)x 10 6 K, probes the bulk of the gas in 
the emitting region. 

• With 400-600 line emission systems detected per 
square degree, one can trace the spatial distribu- 
tion of the line emitting g large fraction of 
which is in the warm-hot phase. The visual in- 
spection of the emission maps confirms that the 
selected line systems trace the large scale filamen- 
tary structure of the gas with typical overdensity 
Pgas/ (j°gas) > 75. Since line emission is associated 
to the denser part of the WHIM, the tracing is 
rather sparse and concentrated around the virial- 
ized structures that punctuate the filaments. More- 
over, because of the unfavorable ionization balance, 
the tracing is biased against the largest virialized 
structures such as galaxy clusters that are typically 
found at the nodes of the network of filaments. 

• The observed dN/dz relation can be compared 
to theoretical predictions to discriminate among 
WHIM models and constrain the physical processes 
that regulate its thermal state and chemical com- 
position. As an example we have shown that the 
dN/dz statistics is capable of discriminating among 
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different metaliicity models with a significance of 
3-5(7 with a 1 deg 2 observation, depending on the 
energy resolution. Further constraints will be pro- 
vided by the analysis of the 3D distribution of the 
line emission systems and by combining these re- 
sults with those obtained from the analysis of the 
absorption spectra. 
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